Canyon of Current Suppression in an interacting two-level Quantum Dot 
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Motivated by the recent discovery of a canyon of conductance suppression in a two-level equal 
spin quantum dot system [Phys. Rev. Lett. 104, 186804 (2010)] the transport through this system 
is studied in detail. At low bias and low temperature a strong current suppression is found around 
the electron-hole symmetry point independent of the couplings, in agreement with previous results. 
By means of a Schrieffer-Wolff transformation we are able to give an intuitive explanation to this 
suppression in the low-energy regime. In the general situation, numerical simulations are carried 
out using quantum rate equations. The simulations allow for the prediction of how the suppression 
is affected by the couplings, the charging energy, the position of the energy levels, the applied bias, 
and the temperature. We find that away from electron-hole symmetry, the parity of the couplings 
is essential for the current suppression. It is also shown how broadening, interference, and a finite 
interaction energy cause a shift of the current minimum away from degeneracy. Finally we see how 
an increased population of the upper level leads to current peaks on each side of the suppression 
line. At sufficiently high bias we discover a coherence-induced population inversion. 

PACS numbers: 73.63.Kv,73.23.Hk,73.23.-b 



I. INTRODUCTION 

The conductance of quantum dots is dominated by 
their discrete level spectra^^ Particularly interesting 
features can be found at degeneracies of the discrete 
levels^— At zero magnetic field spin-degenerate levels 
can be modeled with the Anderson Hamiltonian. 7 At low 
temperatures such systems exhibit the intriguing phe- 
nomenon of the Kondo effect, which has been exten- 
sively investigated^ On the contrary, much less is known 
about systems with orbital degeneracies. The simplest 
system with such degeneracies is the two-level spinless 
quantum dot, which was studied theoretically in connec- 
tion with phenomena such as phase lapses of the trans- 
mission phased— charge oscillationS ) 12 i 13 and correla- 
tion induced conduction resonances^ Relations between 
these phenomena were discussed in Refs. [Ull While 
Refs. I9l4l6l investigated the properties of the system in 
the low-bias limit, the present paper deals with the ef- 
fects of finite bias and temperature. For typical param- 
eters applicable to experiments, we find that changing 
the bias might lead to qualitatively very different results. 
Finite bias was previously investigated in Ref. 1171 where 
the current and full counting statistics for the system 
were studied, but only to lowest order in lead-dot cou- 
pling. The effect of finite bias and temperature in the 
co-tunneling regime has thus not previously been inves- 
tigated. A two-level quantum dot system, where the two 
levels were coupled to separate source and drain con- 
tacts, were studied in a recent paperJ^ Here the orbital 
quantum number has the same effect as spin, due to the 
separate contacts, and orbital Kondo effect can be in- 
vestigated. The paper studied the level renormalization, 
Kondo temperature, local density of states and conduc- 
tance of the system. 

On the experimental side, the degeneracy of orbital 
levels was recently studied for a gate-defined quantum 



dot in an InSb nanowire^ where a canyon of conduc- 
tance suppression was found. Here large level- dependent 
g-factors enabled the study of degenerate orbital levels 
using the Zeeman effect^ Similar systems of quantum 
dots embedded in nanowires have previously been exper- 
imentally realized^ - — The model is not restricted to the 
two-level quantum dot but can also be applied to parallel 
quantum dots coupled to the same leads £2r— 

Previous theoretical investigations of the equal spin 
two-level quantum dot showed the existence of a com- 
plete conductance suppression around the electron-hole 
symmetry point in the limit of zero temperature^^— iSl 
In this paper we investigate the current suppression at 
finite bias and temperature, as well as the current away 
from electron-hole symmetry, i.e. regimes not previously 
studied beyond sequential tunneling. The aim of this 
paper is to determine whether the canyon of current sup- 
pression is generic or appears only for certain parameters. 
In the regime T > Vbias, ksT, where T is the level broad- 
ening given by the coupling between dot levels and leads, 
we find a complete current suppression close to degen- 
eracy, both in the co-tunneling and the sequential tun- 
neling regime, assuming that the two dot states couples 
to the leads with different parity. For weaker couplings 
to the leads, or for increased Vbias or fc^T, only par- 
tial current suppression is found. Results are compared 
with three limiting cases of bias, i) The low-bias limit, 
Vbias < r, U, is shown in Fig. Ufa), where U is the charg- 
ing energy, ii) The high-bias limit F, AE < Vbias < U, 
see Fig. []Jb), where AE is the splitting between the two 
levels, so that both E\ and E% are entirely inside the 
bias window, iii) The ultra-high-bias in which the bias 
is the largest parameter T, AE, U < Vbias; as shown in 

Fig. He). 

For the high- and ultra-high-bias limit the current can 
be evaluated by analytical means, as first-order approxi- 
mations in tunneling are valid. We find that in the high- 
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FIG. 1: Level configuration at low bias (a), high bias (b) 
and ultra- high bias (c). Due to the Coulomb interaction the 
energies of the levels are increased by an amount U when both 
states are filled. The applied bias is given by fiL — (J.R, where 
fiL and fiR axe the chemical potentials in the left and right 
contact, respectively. 



bias limit, the current suppression is completely indepen- 
dent of the couplings . 28 ' 29 In the ultra-high-bias limit, on 
the other hand, only partial suppression is found when 
the two states have different coupling strengths to the 
leads. 

The remainder of the article is organized as follows. 
The system Hamiltonian and the modeling of transport 
is discussed in Sec. II. The canyon of current suppression 
is studied in Sec. Ill using the second order von Neu- 
mann method (2vN), where co-tunneling and coherence 
are included^ The findings of this section are further ex- 
plained in Sections IV- VII. In Sec. IV an intuitive expla- 
nation for the conductance suppression is given by means 
of a Schrieffer- Wolff transformation^ In Sec. V the cur- 
rent at electron-hole symmetry is analyzed as a function 
of bias, temperature, and couplings. The current away 
from electron-hole symmetry is studied in Sec. VI. In 
Sec. VII focus is changed to the current peaks surround- 
ing the canyon of current suppression. At increased bias 
we discover an unexpected population inversion and in- 
vestigate its importance for the peaks. Here the term 
population inversion means a higher population of the 
upper level. It should be pointed out that this popula- 
tion inversion results from a finite bias and is different 
from what is studied in Refs. 12.13. 



II. MODELING OF TRANSPORT 



In this section we discuss the Hamiltonian of the spin- 
less two-level quantum dot as well as the modeling of the 
current through the system. In the single-particle eigen- 
basis of the dot, the system Hamiltonian is given by: 



H — i?dot + -ffleads + Ht, (1) 

H dot = E 1 d\d 1 +E 2 d\d 2 + Ud\d 1 dld 2l (2) 

Pleads = E k c\ p Ckl, (3) 

k,l=L/R 

Ht = (tnd\+t e2 dl)c ke +rl.c, (4) 

k,£=L/R 

where we have assumed that the couplings tu are inde- 
pendent of k and Yn{E) — 2irt^ i po with a constant den- 
sity of states po(E) = J2 k S ( E k - E), for -D < E k < D. 
In the simulations a large bandwidth D is used, assum- 
ing wide conduction bands of the leads. The operators di 
{d\) and c k i (c£j) are annihilation (creation) operators of 
electrons in the dot and leads, respectively. In Eq. @ the 
charging energy U is due to Coulomb interaction between 
the electrons when both dot states are filled. 
The energy levels are parametrized by 



AE 

E 1/2 =±—-E g -U/2. 



(5) 



Electron-hole symmetry, where the levels are symmetri- 
cally placed around the Fermi energy, is given by E g = 
and AE — 0. At this point electrons and holes contribute 
equally to the current through the quantum dot. Bias is 
applied symmetrically to the quantum dot so that the 
chemical potentials are (i^ = Vbias/2 and fiji — — Vbias/2, 
for the left and right contact, respectively. 

In this paper we mainly study couplings of the type 
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where the asymmetry parameter a is chosen to be real. 
Time-reversal symmetry assures that the coupling ele- 
ments are real, implying that the two dot states couple 
to the leads with the same or opposite parity^ This pa- 
per focuses on the case of opposite parity as it is essential 
for the canyon of current suppression at finite bias, see 
Sec. Ill where also the effects of asymmetric couplings 
to left and right leads are briefly discussed. Further- 
more, it was shown in Ref. [l^ that good agreement be- 
tween theory and experiment could be achieved for real 
couplings. As the broadening of the levels is propor- 
tional to the square of the coupling elements we have 

Tli = r^i = r 1: = r^ 2 = o Ti. 

The current through the dot is calculated using the sec- 
ond order von Neumann (2vN) method^ This method is 
an extension of the generalized master equation. Starting 
from the von Neumann equation an equation of motion is 
derived not only for the reduced density matrix but also 
for the elements of the total density matrix consisting 
of a single electron-hole excitation in the contact. In the 
latter system of equations two electron- hole excitation el- 
ements, corresponding to co-tunneling, enter. In the time 
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evolution of the co-tunneling processes there are contri- 
butions from three electron-hole excitation elements but 
these are ignored. This results in a closed set of equa- 
tions, which can be solved for the occupations of the dot 
levels, the coherences between these resulting from the 
coupling to the leads, and the current flowing through 
the quantum dot. 

There are currently several other methods that account 
for tunneling of higher order than sequential processes. 
The most widely used technique is the generalized master 
equation approach which can be derived in many different 
ways including the real-time diagrammatic techniqu o 32 ' 33 
and the Bloch-Redfield approach originally developed in 
Refs. I3414361 Comparisons of different approaches have 
been performed in Refs. I37ll38l . The detailed relation be- 
tween these methods and the 2vN-method will be the 
subject of future investigations, where we show that the 
2vN-method is similar to other methods of 4 th order in 
the couplings t, but contains certain diagrams up to in- 
finite order. For the Anderson model with infinite U it 
has been shown to give the same current as the resonant- 
tunneling approximation 3 — in the real-time diagrammatic 
approach. — This method includes irreducible diagrams 
with an arbitrary number of correlated tunneling pro- 
cesses, but restricted to one electron- hole excitation at 
any given moment. 



III. THE CANYON OF CONDUCTANCE 
SUPPRESSION 

We proceed by investigating whether the canyon of cur- 
rent suppression observed in Ref. [l9| can be found only in 
a certain parameter space or if it is generic for the stud- 
ied two-level system. We plot the normalized current 
J/Vbiasj which equals the conductance G in the low-bias 
limit. As a start, we restrict to couplings given by Eq. © 
and study two different coupling strengths Ti/ksT = 1 
and Ti/ksT = 4. In both cases the relative coupling 
strength of the two levels is given by a — 0.5. The reason 
for using Ti > ksT is that we want to study the cur- 
rent suppression in the co-tunneling regime. The com- 
plete suppression found close to degeneracy in the case 
of Vi/hsT = 4 is also supported by the analytical argu- 
ments of Sec. IV. Furthermore the 2vN-method is non- 
perturbative and contains some classes of diagrams up 
to infinite order, which makes the analysis in the regime 
r > fc^T more reliable. Furthermore, both low bias 
Vbias/fcsT = 1 and high bias Vbias/fcsT = 15 are studied 
to investigate the effects of finite bias. In the left part of 
Fig. [21 the results of the 2vN simulations for U/k B T = 25 
are shown. In order to study the role of the Coulomb in- 
teraction, corresponding results for U — are shown in 
the right part of Fig. [5J Here we applied the method of 
nonequilibrium Green's functions (NEGF)iSi^i providing 
identical results to the numerically much more involved 
2vN method (for U = 0), for the simple system studied 
here. We believe that this holds also for more complex 
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FIG. 2: Conductance calculated with the 2vN method for the 
interacting system, U = 25 ksT to the left, and noninter- 
acting Green's function result to the right. The white areas 
in the right figures are the regions where the noninteracting 
formalism is not valid. The upper four figures correspond to 
the low-bias limit Vbias = ksT, while the lower four show the 
results for high bias Vbias = 15 ksT. The couplings are given 
by Eq. © with Ih = k B T and a = 0.5 for figures (a), (b), (e) 
and (f), while Ih = 4 ksT and a — 0.5 for figures (c), (d), 
(g) and (h). The dashed boxes in the noninteracting figures 
mark the studied region in the interacting case. 
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systems, see also Refs. [301I421 

In order to facilitate the comparison, we performed a 
cut of the noninteracting results at E g — and shifted 
the upper/lower part by ±12.5fcsT. This mimics the 
presence of Coulomb repulsion by one entirely filled or 
completely empty level, as expected in the respective en- 
ergy range. Plotting the noninteracting results in this 
way we find indeed rather good agreement with the in- 
teracting results in the low bias regime, allowing us to 
explain the essential features in this region. For high 
bias the current suppression is more pronounced in the 
interacting case. 

The Coulomb blockade regime present in the left fig- 
ures is a result of the interaction-related correlations. 
This regime can naturally not be studied using non- 
interacting methods. 

The noninteracting results for weak bias, Fig. IHb,d), 
show a line of current suppression, which is bowed to- 
wards the weaker resonance. This can be understood in 
terms of interference between the two levels. The Brcit- 
Wigner formula provides 



t(e) = rf 



f 



E-Ei + iVi E- E 2 + ia 2 T! 
and the current can then be calculated using 



1 



[f L (E) - fn(E)]T(E)dE. 



(7) 



(8) 



Vanishing conductance in the low temperature limit cor- 
responds to T(0) = 0, i.e. complete conductance sup- 
pression can only occur for E2 = a 2 E\. This explains 
the location of the strong current suppression found in 
Fig. [UJb) and (d). It is easily shown that if one level 
couples symmetrically to left and right, complete sup- 
pression is only found for couplings like Eq. ([6]), where 
the sign of tu and tm differs for i = 1,2. This motivates 
the special attention given to these couplings. 

As the bias is increased the current not only depends 
on T(0). When the bias is larger than the width of the 
dip, aT\ at E\ = E 2 = resulting from Eq. ([7]), this 
results in only partial current suppression as can be seen 
in Fig. Eft h). 

Fig. shows that the weak coupling and low-bias 
interacting result is well reproduced by the noninter- 
acting result, Fig. [2jb), in its region of applicability 
(| Eg \> U/2), apart from a slight shift of the suppres- 
sion line. In the Coulomb blockade regime (| E g |< U/2) 
almost no structure is visible due to the very weak co- 
tunneling. In Sec. V we see that this corresponds to the 
regime where ksT dominates and only partial conduc- 
tance suppression can be found. 

Fig. [2fc) shows that the 2vN-method can sometimes 
give unphysical results corresponding to small negative 
currents. This is expected as the 2vN method ne- 
glects some correlated transitions of more than two elec- 
trons/holes. The contributions of these terms become 



important in the limit of large couplings and small tem- 
perature and bias. Apart from this, the 2vN method 
agrees well with the noninteracting results. However, we 
see that the shift of the suppression line is larger when 
the couplings are strong. As explained in Sec. VI this 
can be attributed to the broadening of the levels. In the 
Coulomb blockade regime a clear conductance suppres- 
sion is found. This is the regime where T dominates, i.e., 
the situation described by Refs. [T3 - [l6l 

In the high-bias limit, the interacting and noninteract- 
ing results differ significantly at the canyon of current 
suppression, as seen in Fig. (5Je, g) and Fig. [2Jf, h), re- 
spectively. In the sequential tunneling regime a strong 
current suppression is found in the interacting case. This 
result has been explained by Refs. I28U29I and is further 
discussed in Sec. VI. This interaction-induced suppres- 
sion is not included in the NEGF-method, resulting in 
only partial suppression, as can be seen from the gray 
colored canyon. 

The current suppression in the co-tunneling region be- 
comes weaker as we are approaching the region where 
Vbias dominates, which results in only partial current 
suppression, see Sec. V. This shows that the canyon of 
current suppression is found in all four cases albeit the 
current not always entirely drops to zero. Further in- 
vestigations show that the canyon of suppression is not 
restricted to couplings like Eq. but is present also 
for couplings which are asymmetric with respect to the 
left and right contact. The importance of the relative 
coupling strength, a, is investigated in Sec. V. 
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FIG. 3: Conductance calculated with the 2vN method. The 
two levels have equal parity with couplings given by tm = 
tm = t, tL2 = 0.4t and tR2 = 0.6i with Tli = r^i = 
2npot 2 = 5ksT. Additional parameters are Vbias = 15 fesT, 
U = 25 k B T. 

For couplings where the two levels have the same par- 
ity, i.e. the two levels couple with the same sign to 
left and right lead, the situation is different, as shown 
in Fig. [3] Here we still see the high-bias blockade in 
the sequential tunneling regime i 28 ' 29 For E g < —U/2 or 
Eg > U/2 interactions are of less importance and the 
blockade disappears. One also sees that no canyon is 
visible in the Coulomb blockade regime. The reason is 
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that the width of the dip in the transmission function be- 
comes increasingly narrow as the singular coupling point 
is approached. At such couplings, one can at degeneracy 
completely decouple one level from the leads, resulting 
in a bistable system, see also Refs. [I3fl5l As a result, 
the canyon is present only for very low bias and temper- 
ature, i.e. a regime where the 2vN-method is not appli- 
cable. Fig.|3]clearly shows that at finite bias the parity of 
the couplings is essential for the suppression, unlike the 
zero-bias limit. Two levels with the same parity will not 
be studied any further. The remainder of the paper is 
devoted to the couplings described by Eq. ©. For such 
couplings, details of the canyon of current suppression is 
summarized in Table I, which can be of use for father 
experimental studies. 

TABLE I: A summary of the current suppression from Fig. [2] 
For levels with couplings of different parity the current sup- 
pression is generic, but not always complete. The high-bias 
limit of Fig. \V[b) implies to E g ~ ±(7 /2. Therefore the co- 
tunneling regime is not applicable (n.a.) in this limit. In the 
ultra-high-bias limit the current is given by the sequential 
tunneling result. The co-tunneling regime is thus not catego- 
rized in this limit. The table also indicates in which sections 
the different regimes are investigated further. 



r 


Vbias 


co-tunneling 
(E g » 0) 


sequential tunneling 
(E g « ±U/2) 


> k B T 


low 


complete 
(Sec. IV and V) 


complete (Sec. VI) 


> k B T 


high 


n.a. 


complete (Sec. VI) 


> k B T 


ultra-high 




partial (Sec. V) 


< k B T 


low 


partial (Sec. V) 


partial (Sec. VI) 


< k B T 


high 


n.a. 


complete (Sec. VI) 


< k B T 


ultra-high 




partial (Sec. V) 



the Hamiltonian in terms of pseudo-spin operator a 15 i 16 
defined by = dx, d± = d 2) c k \ = (c kL + c w )l^ and 
c/4 = {—CkL + Ckn)/\/2. The system Hamiltonian then 
reads 

H dot = ^ t 4d t + ^4^ + C/ 4 d t4 d i' ( 9 ) 

Pleads = E k C L C k°, (10) 

H T = ^X4 CT <2 ff +H.c. 

ak 

= V^(4t d t + ac^dx) + H.c. (11) 
k 

Here E t = Ex, £4. = E 2 , t t = \/2t, t± = V%at, and the 
problem has been mapped onto the one-lead Anderson 
model with spin-dependent couplings. A similar trans- 
formation can be found for arbitrary couplings. 

We now perform a Schrieffer- Wolff transformation, i.e. 
a canonical transformation 

H s = e lS He~ lS , (12) 

so that Hs does not have a linear term in ■ One finds 
that the following transformation S = S + + S~ , with 

+^ 2 ^-(i-^)4^) I as) 

and S + = (S'~) t , will do the job. Here ~5 =\. if a =f 
and (7 = j" if a Writing the tunneling Hamiltonian as 
H T = H+ + H T , with 



IV. TUNNELING PROCESSES IN THE 
LOW-ENERGY LIMIT 

In this section we study the current in the regime where 
E x , E 2 < and Ex + U, E 2 + U > 0, so that the 
dot is singly occupied. Furthermore we will assume that 
Vbias, ksT — > 0. This allows us to derive an effective low- 
energy Hamiltonian by means of a Schrieffer- Wolff trans- 
formation. Low-energy implies that all tunneling events 
occur around Eu — 0. The transformation requires that 
the charging energy of the dot is much larger than the 
broadening of the dot levels, so that the dot is singly 
occupied when the levels are placed as in Fig. [IJa) . By 
analyzing the effective Hamiltonian in this regime we will 
be able to identify the two-particle processes contribut- 
ing to the current and give an intuitive explanation of 
the conductance suppression. It is convenient to write 



= Y t tA* d *i ( 14 ) 

ak 

and — [H^,)\ the effective low-energy Hamiltonian 
of the singly occupied subspace, when restricting to two- 
particle tunneling, is given by 

Hs = -ffdot + -Hleads 

+ % -{[S-,H+] + [S + ,Ht))+0{H*) (15) 

Hts 

where the last term is the effective tunneling Hamiltonian 
Hts- Performing the commutations and returning to our 
original basis one finds: 
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1 ( t 2 dU 2 t 2 il- dU 2 )\ , <r + \ / 

H ™ = 2^ \ E k -(E l+ U) + E k E[ ) + C ^ )iCk ' L + Ck ' R) 

+ 2E U-(£b + t/) + i-E, J (C - " C - )(C - " (16) 



2 L + " J (Cfei + cL)(Cfc,£ " Ck,R)d2dl 

kk' v v ' / 



+ H.c. + F, 



where indicates that the terms containing ct L CkL + 
c\ R CkR, corresponding to the renormalization of energy 
levels -ffrcnorm, are not included in the sum. In the limit 
of zero temperature and zero bias, the renormalization 
reads (see Appendix B) 



fr — El I In 

^rcnorm I rn 

7T 



-a 2 ln 



Eh 



Ei + U 
E 2 



d\di 



E 2 + U 



d\d 2 



(17) 



It should be noted that the above expression does not 
hold if one of the levels is close to the chemical poten- 
tial, as the Schrieffer- Wolff transform breaks down in this 
regime. The effect of the renormalization will be further 
discussed in Sec. V, here we simply note that at electron- 
hole symmetry, where E\ — E 2 — —U/2, the renormal- 
ization vanishes. 

Each term in Eq. (|16j) corresponds to a co-tunneling 
process, where the electron is transferred from one lead 
to the other, or returns to the original lead. The first two 



sums correspond to elastic processes where the state of 
the dot is unchanged. The next two sums correspond to 
inelastic processes where the state of the dot changes. 

Generally there is an effective overlap matrix element 
between the two dot states originating from the coupling 
to the leads. In Refs. [27ll43l it was shown that for cou- 
plings of the type Eq. ([B]), this overlap vanishes. Indeed, 
evaluating the k = k' inelastic terms in Eq. (1161) results 
in factors c kL CkL — c \u c kR = in the low-bias limit, see 
Sec. VII for a discussion concerning finite bias. 



Let us now analyze Eq. (|16|) in detail: We see that 
only the inelastic processes affect the occupations of the 
dot, as d\d\ and d\d 2 commutes with the other parts of 
Hts- For E\ = E 2 the third and fourth sum in Eq. (TT6"1) . 
corresponding to the inelastic processes, are symmetric in 
level 1 and 2 except for a sign change ^fc^Ck' n±ct R Ck> l- 
In the low- bias limit, this sign change does not affect the 
kinetics of the occupations and we conclude that both 
levels are half filled at degeneracy. 

The current operator reads Js = — Ylk tl^s> c ! c fc]> 
resulting in 



Js 




d\d 2 



+ 



d\di 

{E 1 + U) ' Ek-Et 
1 1 



E k - {E 1 +U) E k - E x 



a 2 d\di 



i 2 dtd? 



Ek — (E 2 
d\d\ - 



U) 



2 2 \ / t t \ 

Ek ^~E~2 J (c ^ Cfc ' L ~ c kL c k'R) 



(18) 



H.c, 



Ek - {E 2 + U) Ek- E 2 



(cl R Ck>L + c\ L Ck'R) 



where d\d\ + d\d 2 — 1 has again been used. The first 
row corresponds to elastic processes, while the second 
row constitutes the inelastic processes. We adopt the no- 
tation | 0), | 1), | 2), and | d), corresponding to the dot 
being empty, in state 1, in state 2, or doubly occupied. 



The term d\d 2 j(Ek — {E\ + U)) represents the elastic pro- 
cesses | 2) — s- | d) —> | 2). The denominator includes the 
energy difference between the two states involved in the 
process, E d - E 2 = E x + E 2 + U - E 2 = E x + U, and the 
numerator represents the fact that state | 2) has to be oc- 
cupied for this process to occur. In the same way the term 
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d\di/(Ek — Ei) represents the process | 1) — > | 0) — > | 1). 
Assuming Ek — 0, corresponding to low bias and low 
temperature, we see that these terms cancel completely 
at the electron-hole symmetry point E\ = Ei = — U/2, 
owing to the equal population of the two levels. We note 
that these two processes, involving tunneling into and out 
of state | 1), use the couplings of this state. In the same 
way we see how the two processes | 1) — >■ | d) — > \ 1) and 
| 2) —> | 0) —> | 2) cancel. The factor of a 2 enters here as 
it is state | 2) that is active in these processes. We have 
seen that the elastic co-tunneling processes cancel com- 
pletely at the electron- hole symmetry point, independent 
of the relative coupling strength of the two levels. Away 
from electron-hole symmetry, we see that the processes 
no longer cancel. However, at E\ = E%, we again find 
complete canceling for a — I. 

We now investigate the inelastic processes. It can eas- 
ily be seen that the terms involving d 2 d\ and d\d2 cancel 
at Ei = E 2 , due to the equal occupation of the two lev- 
els. Electron- hole symmetry is thus not required for the 
canceling of the inelastic terms, degeneracy of the two 
levels is sufficient. 

To summarize we have shown that at Ei = E% = —U/2 
the current suppression is complete, in the limit of low 
bias and low temperature at the electron-hole symmetric 
point. Unlike previous explanations of this phenomena 
our derivation gives an intuitive explanation to the sup- 
pression in terms of canceling of co-tunneling processes. 
It also explains why the suppression is not complete for 
finite bias and temperature or away from electron-hole 
symmetry. This will be further studied in Sees. V and VI 
respectively. The suppression is not an interference-effect 
between different processes, since the processes that can- 
cel involve different initial and final states. The pro- 
cesses instead cancel as a result of the equal population 
of the two levels. At degeneracy, but away from electron- 
hole symmetry, the elastic processes do not cancel unless 
a = 1. There is, however, partial canceling due to the 
different parity of the two levels. We stress that the 
complete suppression at electron-hole symmetry is not 



an effect of this difference in parity. The pseudo-spin 
operators are always orthogonal. It is this orthogonality 
that gives equal occupations of the two levels and ensures 
the canceling of current carrying processes in the low-bias 
limit. The canceling is due to the phase lapse between 
the elastic co-tunneling processes using | 0) and | d) as 
the virtual intermediate state. This shows the relation- 
ship between population switches, phase lapses, and van- 
ishing conductance, which was previously investigated in 
Ref. lla. At electron-hole symmetry the vanishing con- 
ductance is thus an effect of the correlations induced by 
the Coulomb interaction. 

V. CURRENT SUPPRESSION AT 
ELECTRON-HOLE SYMMETRY 

In the previous section we saw that the current sup- 
pression was most pronounced at the electron-hole sym- 
metry point. In this section we will investigate the sup- 
pression at this point in more detail. The current is cal- 
culated as the flow of particles from the left contact into 
the dot. Due to current conservation this is, in the sta- 
tionary case, of course equal to the current flowing from 
the dot into the right contact. Thus, we do not label the 
current with any lead-index. To quantify the extent of 
the current suppression we introduce the dimensionless 
ratio between the current through the two-level dot cor- 
responding to Eq. ([2]) and the current through a single 
level dot containing only level one 

Q^both levels /-i n\ 

= j ■ l iy J 

<>levcl 1 

This ratio measures how much the presence of level 2 sup- 
presses the current. In the regime where sequential tun- 
neling is dominating the first order von Neumann method 
(lvN) can be used to approximately evaluate the current, 
see Appendix A. Using the couplings of Eq. ([6]), it is 
possible to derive Eq. (|A8[> . At electron- hole symmetry, 
Ei = E 2 = —U/2, this simplifies to 



J = {w 00 f L (-U/2)(l + a 2 ) - w n [l - fd-U/2)] +wnf L (U/2)a 2 - w 22 [l - f L (-U/2)]a 2 
n 

+W22MU/2) - w dd [l - f L (U/2)](l + a 2 ) + 2aTZ{wi 2 } [1 - f L (-U/2) + f L (U/2)}} , (20) 

I 



where woo, u>n, W22 and Wdd is the probability to find 
the dot empty, in state 1, in state 2 and doubly occupied 
respectively, while W12 is the coherence between level 1 
and 2. Here 1Z denotes the real part. In this paper we 
refer to the coherence as the non-diagonal elements of the 
reduced density matrix. The coherence is thus a basis- 
dependent quantity. The first terms of Eq. (|20p are iden- 
tical to the corresponding Pauli master equation and are 



easily interpreted in terms of transition rates. The last 
term originating from the coherence of the two levels is 
more difficult to understand intuitively. However, it is 
evident that a negative real part of the coherence results 
in a decreased current. 

At level degeneracy one can get rid of the coherence 
term by performing a basis change that diagonalizes the 
reduced density matrix. Here we focus at the case of 
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electron- hole symmetry, E\ = E2 = —U/2, where level 
1 and 2 have the same occupation and the coherence is 
real (see Appendix A). The singly occupied part of the 
reduced density matrix can thus be written 



Wn W12 

w 2 i w 2 2 



wn W12 
w 12 wn 



SDS- 



the applied bias. The numerator of Q is calculated us- 
ing the 2vN method, while the denominator is calculated 
using NEGF which, however, agrees exactly with 2vN 
for the single level dot. Fig. 2] clearly shows that when 
Vbias > U, the ultra-high-bias regime, where sequential 
tunneling is dominating, is entered and the results of the 
2vN simulations agree very well with the analytical re- 
sults, Eq. pgj). 



The diagonalization is given by 




D = 



w\\ + W12 
wn - W12 



where the diagonal elements of D are the occupations in 
the new basis. 

The couplings in the new eigenbasis are given by 



*lx = ^( 1_a )' *Sa = 7f (1 + a) ' 



(21) 



In the new basis one level couples strongly to the left, 
while the other couples strongly to the right lead. This 
results in a strong current suppression and also explains 
the occupations in our new basis: the level which is 
strongly coupled to the lead with a high chemical po- 
tential has a high occupation. A high bias thus results in 
a large coherence between the levels, cf. the expression 
for the matrix D. For a = 1, i.e. both levels have the 
same coupling strength, each level decouples completely 
from one lead. As a result no current can flow through 
the quantum dot, see also Ref. H^. The effect of the 
coherence is thus to decrease the current. As the coher- 
ence is zero in the new basis one can use the Pauli master 
equation approach to derive the current. At electron-hole 
symmetry we obtain 



Jb 



oth levels 



(i-a 2 ) 2 r! 



[f L (-U/2)-f R (-U/2)}.(22) 



2h(l + a 2 ) 

If only level 1 was present the current would be given by 



level 1 



£1 
2h 



[f L (-U/2) - f R (-U/2)], 



and the resulting ratio is 



Q = 



■/both levels 
■he 



cl 1 



1 



(23) 



(24) 



This result is valid when sequential tunneling is dominat- 
ing. 

In Fig. 2] the current suppression, Q, at electron-hole 
symmetry, E\ = E2 = —U/2, is shown as a function of 
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FIG. 4: Q as a function of H>ias at electron-hole symmetry for 
different coupling strengths. All curves correspond to a = 0.5 
and U = 25 k B T. 

In the low-bias limit the results are more complicated. 
For weak couplings the contribution from co-tunneling 
processes is small and the first-order result Eq. (|24|) is 
approached. As the bias is increased the phase space for 
co-tunneling processes (~ Vbias) increases, so these can 
no longer be neglected. This results in a stronger current 
suppression and a related decrease in Q for the weaker 
couplings. For strong couplings, i.e. T > Vbias ,kBT, 
co-tunneling dominate and the low-energy Hamiltonian 
of Sec. IV predicts complete suppression in the low-bias 
limit. We stress that the suppression in the low-bias limit 
is not an artifact of neglecting higher order processes. 
As was shown in Sec. IV and Refs. QJ-QJ and |H the 
suppression is complete in the limit of low bias and low 
temperature, at the electron-hole symmetric point. 

In Fig.[5]the simulated current suppression at electron- 
hole symmetry is shown for different temperatures as a 
function of bias. As the temperature is not constant in 
Figs. 5 and 6 we here express all energies in terms of Ti. 
It is clearly seen that the sequential tunneling result is 
reached in the high-temperature limit. 

In Fig. [6] we investigate the current suppression Q as a 
function of the difference a in coupling strength between 
the two levels for different temperatures. For lower tem- 
peratures it can be seen that there is essentially com- 
plete suppression when both levels have equal coupling 
strength, i.e. awl. When a w and the difference in 
coupling strength is large, almost no suppression is ob- 
served. This is expected as no current flows through the 
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FIG. 5: Q as a function of Vbiaa at electron-hole symmetry 
for different temperatures. All curves correspond to a — 0.5 
and U = 5 IV 
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FIG. 6: Q as a function of the difference in coupling strength, 
a see Eq. (J6j , at electron-hole symmetry for different temper- 
atures. The high temperature limit is given by Eq. (|24|) . All 
curves correspond to Vbias = 1 fi and U = 5 IV 



weakly coupled state, and due to electron-hole symmetry 
it does not matter if it is occupied or not. 

In this Section we investigated how a finite bias and 
temperature affects the suppression at electron-hole sym- 
metry. In the limit of low/high bias and temperature we 
compared with previously established results/first order 
simulations, and found good agreement. In intermediate 
regimes we found that the current suppression could show 
a non-monotonous behavior, see Ti = 1.25 ksT in Fig. 4. 
Furthermore Figs. IH [5] and |6] show that Q increases with 
increasing temperature until the high temperature limit 
is reached. We have seen how the current suppression is 
given by analytical first order results in the regimes where 
sequential tunneling is dominating. In this regime the 
suppression is coherence-induced and only partial unless 



both levels have the same coupling strength. Stronger 
couplings lead to an increased co-tunneling which results 
in a stronger suppression. 



VI. CURRENT SUPPRESSION AWAY FROM 
ELECTRON-HOLE SYMMETRY 

With a clear picture of the current suppression at 
electron-hole symmetry, we now change focus to the 
regime away from this symmetry point. We will inves- 
tigate the current at E g = — {7/2. At this point the 
singly occupied levels are close to chemical potentials of 
the leads, and sequential tunneling contributes strongly 
to the transport. Thus, the lvN-method gives reason- 
able results and we can compare first- and second-order 
von Neumann to investigate the importance of the broad- 
ening in this regime. Fig. E£a) shows the current at 
Eg = -U/2 = -12.5 k B T calculated with 2vN and lvN 
methods. Since level broadening is not included in the 
first order method the peaks are too high and narrow. 

More interestingly, we see that the current minimum 
has been shifted away from degeneracy. Comparing with 
the noninteracting results of Fig. [21 where the dip is found 
at AE — at E g — ±U/2, one draws the conclusion that 
the shift is related to the finite charging energy U. We 
will argue that this shift has three origins: 

i) A small part of the shift originates from a first order 
level renormalization and a bias induced anticrossing be- 
tween the levels, as can be seen from the slightly different 
dip location of the two lvN-curves in Fig. [7^,. The first 
order simulations not including these effects also show a 
shift, however slightly smaller. This suggests that there 
are additional mechanisms producing the shift. The bias 
induced anticrossing is further discussed in Sec. VII in 
connection with the population inversion. 

ii) The main shift within the first-order approach can 
be attributed to the asymmetry in couplings, which has 
to be compensated by different spectral weights in order 
to allow for full canceling of the individual channels. As 
discussed in Ref. [l2| the spectral density of each state is 
divided among two peaks, situated at E and E + U . In- 
creasing the occupation of one level results in a shift in 
the spectral function of the other level towards the up- 
per peak at E + U. At E g = —U/2 the upper peak is 
far from the Fermi level of the leads and does not con- 
tribute to the transport through the dot. In order for the 
stronger/ weaker coupled level (i.e., levels 1/2 for \a\ < 1) 
to have lower/higher spectral weight at the Fermi en- 
ergy, the weaker coupled level must have a higher occu- 
pation. This is naturally the case, if this level is lower 
in energy, thus at E g w —U/2 the current minimum re- 
quires a positive AE, c.f. Eq. ((SJ. This effect is re- 
duced for increased bias. For Jl{Ei) = fL{E%) — 1 and 
fn(Ei) = fR(E 2 ) — 0, complete conductance suppres- 
sion can be found at degeneracy due to an increased oc- 
cupation of the weakly coupled level, as discussed below. 

iii) Comparing the first- and second-order methods it 




FIG. 7: (a) The current through the two level quantum dot 
calculated with 2vN and lvN formalisms, with and with- 
out renormalization of the energy levels, at Ti = 4 k B T 
U = 25 k B T, E g = -(7/2, a = 0.5 and V hias = 2.5 k B T. As 
the broadening is not included the first order method over- 
estimates the peak height but underestimates the shift away 
from degeneracy of the conductance dip. (b) and (c) show the 
corresponding occupations and and real part of the coherence 
calculated with the 2vN formalism. The vertical dashed line 
marks the location of the conductance minimum. 



is clear that the inclusion of broadening of the levels in- 
creases the shift of the conductance dip. It is well known 
that in a broadened level the low energy part of the spec- 
tral density is mostly occupied^ Due to the Coulomb 
interaction this results in a significantly higher occupa- 
tion of the strongly coupled level at level degeneracy, see 
Fig. EJb) and Refs. [l3lfl3 . To compensate, the minimum 
is shifted so that the strongly coupled level has a higher 
energy than the weakly coupled. Due to the level broad- 
ening one can not expect the two levels to contribute 
equally to the current away from the degeneracy point. 
As the levels are equally far from the bias window, the 
broader level will have a stronger contribution if the oc- 
cupations of the two levels are equal. Thus the minimum 
can be found for values of AE such that the occupation 
of the weakly coupled level is slightly higher, in agree- 
ment with Fig.[7]Jb). In the regime | E g |> U/2 Coulomb 
interaction is of minor importance but the dip is still 
found where the occupation of the weakly coupled level 
is higher, to compensate for the larger broadening of the 
strongly coupled level. This implies the strong shift out- 
side the Coulomb-blockade region, as obtained by NEGF. 

Note that the above explanation holds not only in the 
sequential but also in the co-tunneling regime — U/2 < 



Eg < U/2, but the argumentation has to be reversed for 
Eg > 0, where the current is mainly carried by holes 
through the doubly occupied states. At electron-hole 
symmetry the singly and doubly occupied levels have the 
same distance from the bias window, which results in 
minimal conductance at level degeneracy. 

The canceling of the different current paths through 
the dot is related to coherence between the levels. For 
the couplings considered here we concluded in Sec. V that 
a negative real part of the coherence led to a decreased 
current. It is thus not surprising that Fig. [TJc) shows 
that the current minimum coincides with a minimum in 
the real part of the coherence. 

The relation between the charging energy and the shift 
of the dip can be better understood by studying the cur- 
rent at Eg = —U/2 for different values of U. This is done 
in Fig.|8]for Vbias = 2.5 fc^T. Here U = corresponds to 
electron-hole symmetry and the conductance suppression 
is not complete due to the finite bias. As the charging 
energy increases, the current will mainly flow through 
the singly occupied dot states. As predicted above, this 
results in a shift of the conductance dip towards lower 
energy of the weakly coupled state. Due to the different 
occupations the current through the two levels cancel for 
larger U, resulting in a strong conductance suppression. 
Thus a larger U results in stronger conductance suppres- 
sion at the dip, in agreement with Fig. |HJ 

We now focus on the case of higher bias. When both 
levels are in the bias window and fcsT <C Vbias the current 
can be calculated using the first-order results. Here it 
is important whether the charging energy is dominating 
over bias or not. For U -C Vbias, i-e. in the ultra-high- 
bias limit [Fig.OJc)], both the singly and doubly occupied 
levels are inside the bias window and the current can be 
calculated using Eq. (IA17[) . 
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FIG. 9: The conductance of the two level quantum dot at 
ri = 5 k B T, a = 0.5, V bias = 30k B T and E g = -(7/2. As 
predicted, increasing U leads to reduced conductance at level 
degeneracy. 



were close to the bias window, sequential tunneling was 
important, and different parity of the two levels was es- 
sential for the suppression. At finite bias, the suppression 
was not complete unless the levels had the same coupling 
strength. However, a finite Coulomb interaction shifted 
the dip away from degeneracy and resulted in a stronger 
suppression. This shift, not previously discussed in the 
literature, was studied in detail. Three origins for the 
shift were found: (i) a renormalization of the energy lev- 
els, (ii) a shift required to give the weakly coupled level 
a higher occupation which reduces the spectral weight of 
the singly occupied state for the other level in the case of 
Coulomb repulsion, and (iii) a shift related to the broad- 
ening of the levels. Furthermore, in the high-bias regime, 
Fig.[TJb), there was complete current suppression at the 
degeneracy point. 



VII. CURRENT PEAKS AND POPULATION 
INVERSION 



(l + a 2 )[(l-a 2 f + (Ei-E 2 f/rj] 
h 2[(1 + a 2 ) 2 + (E 1 - E 2 y /T\] ■ 1 ' 

in agreement with Ref. \XA The dip is found at E\ = 
E 2 but the suppression is not complete. (Note that this 
result corresponds to electron- hole symmetry.) 

For U — > oo, i.e. the high-bias case [Fig. Hfb)], the 
doubly occupied state will be empty and the current is 
given by Eq. (|A22|1 . 

ri (i + a^m-E^/ri 

h (i + a 2 ) 2 + 3(£! -£ 2 ) 2 /r 2 ' 1 ' 

again in agreement with Ref. [l7|. We note that the 
current is zero at E\ = E 2 . This can be explained 
intuitively^^ At degeneracy we are free to work in 
any linear combination of our original single particle ba- 
sis. Especially we can choose a basis such that one level, 
^filled = afi + (3ty 2 , is completely decoupled from the 
lead with the low chemical potential. This level will al- 
ways be occupied due to the high bias. However, the 
other level cannot be occupied as U — > oo. As a result 
neither level can carry any current through the dot, as 
seen for U = 50 fc^T in Fig.|9l It is evident that if level 1 
is more strongly coupled than level 2, then /3 > a imply- 
ing that the weakly coupled level has a higher occupancy, 
as discussed above. 

At high bias, as in Fig. [9j the conductance minimum 
can be found close to level degeneracy in agreement with 
Eqs. (I2~5j) and (f2l)|). The symmetry with respect to mir- 
roring in AE = 0, for U = is broken by a finite U. 
Comparing Figs. [8] and [9] shows that increasing the bias 
increases the width of the dip and broadens the peaks. 

In this section the current suppression away from 
electron-hole symmetry was investigated. Since the levels 



In this section we change the focus from the dip to 
the surrounding current peaks. In the limit of zero bias 
and zero temperature and wqq — wad = 0, it was shown 
in Ref. [15] that, for the couplings of Eq. ©, the peaks 
correspond to a difference in occupation of \wn — w 22 \ = 
0.5. To verify if the 2vN method can reproduce these 
results we must go to the regime T ^> Vbi aS j ^bT, i.e. to 
the boundary of the region of validity of the 2vN-method. 
To reduce the possible occurrence of negative currents we 
use couplings with a = 1, where the current is inherently 
zero at degeneracy. The result of the simulations is shown 
in Fig.[TUl where the black curve corresponds to low bias, 
Vbias = 2 fcsT 1 . Closer examination reveals that the peaks 
are located at AE = ±3.4 ksT where \wu — w 22 \ s» 0.4 
in good agreement with the analytical prediction. It is 
not surprising that the difference in occupation is smaller 
than predicted since there is also a finite probability that 
the dot is empty or doubly occupied. 

A more surprising result is found when the occupations 
for higher values of Vbias are studied, see Fig. [TT] This 
is a regime of validity for the 2vN method due to the 
increased bias. It should be pointed out that the regime 
described here, where the upper state is far below the 
chemical potentials, can not be studied in pure second- 
order perturbation in T. Such methods do not include 
mechanisms for emptying the upper state d 5 ' 46 On each 
side of the level degeneracy there is a region of popula- 
tion inversion, i.e. the upper level has a higher occupa- 
tion. This population inversion has not previously been 
investigated for our system, however, a similar heating 
of single electron transistors due to inelastic co-tunneling 
processes was discussed in Ref. H3- The 2vN-simulation 
for the single-level dot with spin, where a magnetic field 
splits both spin directions, (Anderson model, not shown) 
reveals no population inversion between the two spin lev- 
els. This suggests that the coherence between the two or- 
bital levels is essential for the inversion. The coherence 
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FIG. 10: Conductance at £" 9 = of the two level quantum 
dot as a function of the level splitting for different values 
of Vbias- The insert shows the occupation of the two single 
particle levels for Vbias = 2 fcsT, with the vertical dashed 
lines marking the location of the peaks. The parameters are 
given by a = 1, Ti = 20 k B T and U = 100 k B T. 
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FIG. 11: Occupation of the two single particle levels for differ- 
ent values of Vbias- Around degeneracy a region of population 
inversion can be found. The parameters are given by a = 1, 
Ti = 20 ksT and U = 100 ksT. The population inversion 
in the diagonal basis (| a), | /3)) is calculated using fi from 
Eq. p7|. 



simply means that an electron can be in a superposition 
of the two states, which results from the overlap between 
the dot levels introduced by the couplings to the leads. 
For the level configuration of Fig. QT] the dot can as a 
good approximation be considered to be singly occupied, 
which means that the Schrieffcr- Wolff transformation of 
Sec. IV can be used to evaluate the overlap. As noted in 
Sec. IV, in connection with Eq. (|16p . the contribution to 
the overlap from state in the leads is proportional to 



bias results in an effective overlap between the two states. 
Replacing the sums with integrals and using the notation 
of Eq. ([5]) one arrives at the following expression for the 
anticrossing ft 



ft = -at 2 
2 



In 



V bias -U-AE 



-fin 



Vbias + U + AE 

Vbi as -U + AE 



U-AE 



(27) 



in the additional term of the dot Hamiltonian given by 
il(d\d 2 + d\di). We have assumed k B T = 0. As long as 
AE < U, which always is the case if Eq. (|2~T1) is valid, 
f2 is negative for Vbi as > 0. For small bias Vbi as < AE 
the anticrossing vanishes. This is expected as inelastic 
processes where the state of the dot is changed is only 
possible as long as Vbias > AE. To investigate the effect 
of this anticrossing we diagonalize the dot Hamiltonian. 
The resulting eigenvalues and eigenvectors are given by 



E a /p ^1(e 1 +E 2T \/AE 2 + 4ft 2 ) 



| a) = x al | 1) + x a2 | 2) , 
| p) = xpx | 1) + x [j2 | 2) , 

where the coefficients are given by 



(28) 



(29) 



11 = 2iT \ AE ~ ^ AE2 + An J ' Xa2 = Na (30) 



2ft 



[AE + VAE 2 +AQ 2 ) , x a2 = N p 



respectively, with A^ Q and being normalization con- 
stants. Clearly we always have E a < Ep. The new eigen- 
states also get new effective couplings 



tha = (x a lt - X a2 at), t Ra = (Xa.lt + X a2 Clt), 

thp = (xpit - x 02 at), t R p = (xpit + xp 2 at). (31) 

Due to ft < both coefficients x a \ and x a2 of the lower 
eigenstate | a) have the same sign. Combined with the 
difference in parity of our two original states this results 
in a weak coupling of | a) to the left lead, i.e. the lead 
with high chemical potential, and a strong coupling to 
the right lead. In the same way | j3) couples strongly to 
the left and weakly to the right, see Fig. [12] The popula- 
tion inversion is real in the sense that it is present in the 
diagonal basis of the dot (| a), \ /?)), see Fig.QTJ It might 
at first seem surprising that the population inversion is so 
large. For V Mas = 10 k B T in Fig. [TT]we have ft « 4 k B T 



c kL c kL — c kR c kR =^ fhiEk) — fn(Ek), which vanishes over the entire range of AE. The detuning Ep — E a = 2ft 

at Ei — E 2 = 0, i.e. it is almost the size of Vbi asi and 
yet a strong population inversion is observed. The cause 



in the low bias limit. As a result the coherence and the 
overlap between the dot state vanishes. However, a finite 
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FIG. 12: The finite bias induces an anticrossing Q between 
the two original states. Diagonalizing the dot Hamiltonian 
results in the upper state being strongly coupled to the lead 
with high chemical potential, while the lower state is strongly 
coupled to the lead with low chemical potential. 
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for this population inversion is inelastic processes driving 
a current from left to right by filling the state | /3) and 
emptying | a) . Due to the asymmetric couplings of the 
(| a), | /3))-basis such processes are very strong. The pop- 
ulations of Fig. [IT] in the (| a), \ /3))-basis are calculated 
assuming f2 given by Eq. (|27[) . Deviations from single 
occupation of the quantum dot, due to the broadening 
of the dot levels and the finite bias, cause slightly dif- 
ferent values of f2. However, at AE — 0, | a) and | /3) 
are the symmetric and antisymmetric combination of | 1) 
and | 2) respectively, independently of fi. Therefore, the 
population inversion is independent of f2 at AE = 0. For 
higher bias this results in populations close to w aa = 



and 



1. In the limit Vyj n , — > oo, it can not be 



considered as a population inversion since il vanishes in 
this limit, resulting in E a = Ep. 

When AE < 0, i.e. E\ < E2, \a) is mainly composed 
of the original state | 1) , while | 0) is mainly composed 
of I 2). For AE > the situation is the reverse. Thus, it 
is always the case that the upper state | f$) , that couples 
strongly to the left, is mainly composed of the upper 
state in the original basis. This explains the observed 
population inversion in our original basis. A negative 
bias results in ft > 0, which leads to a strong coupling 
of I a) to the left. However, for a negative bias the left 
lead has a low chemical potential. Thus the explanation 
holds also for negative bias. 

The k = kl terms of Eq. (|16[) are present also in first 
order perturbation in dot-lead couplings. We have seen 
that the real part of the sum, (or integral in the continu- 
ous limit), over k — k' gives rise to the population inver- 
sion. We thus expect to see population inversion also in 
first order simulations if the real part is included. How- 
ever, first order perturbation theory can not be applied 
to the co-tunneling regime where Eq. (| 16[) is valid. We 
instead investigate the ultra-high bias regime where all 
states are inside the bias window. Eq. (IA2[) can be solved 
including the real part of the sums. The result is shown 
in Fig. [TBJ where also the result of a 2vN-simulation is 
shown as comparison. One observes that as the bias is 
increased, which means that the real parts of the sums be- 
comes smaller, the occupations approach the infinite bias 
result of Eq. (|A17[) . It is evident that the broadening re- 
sults in a smaller population inversion. We can therefore 



FIG. 13: Populations calculated with the lvN- and 2vN- 
methods. Couplings and charging energy are given by Ti = 
5 k B T and U = 100 k B T, respectively. 



not completely rule out the possibility that the popula- 
tion inversion is an effect of neglecting higher order cou- 
plings. As the inversion is present also at very high bias, 
where the contribution from higher order terms is small, 
we believe, however, that it is a real effect. Since it is the 
real part of of the sums of Eq. (IA2[) that is responsible 
for the renormalization one would not expect to see any 
population inversion when these are neglected. Indeed, 
already at Vu as = 200 ksT one finds that the results of 
such simulations are indistinguishable from Eq. (|A17[) . 
These results are therefore not shown in Fig. [TBI 

The population inversion is in fact related to the con- 
ductance peaks, by increasing the elastic co-tunneling. 
Let us assume that E% > E± so that W22 > wn- Then in 
Eq. ([T5| the elastic processes with d\d,2 will dominate. If 
the levels are placed as in Fig. [IJa), then Ek — (E\ + U) 
and —{Ek — £2) are both positive and there is construc- 
tive interference between these tunneling processes. Put 
simply, the population inversion activates elastic pro- 
cesses through the states close to the bias window, i.e. 
processes that contribute strongly to the current. Fur- 
thermore, the elastic processes have the same phase and 
interfere constructively. For large detuning, AE, the co- 
herence decreases and as a result the occupation of the 
upper level will decrease. Increased temperature, which 
also results in a lower coherence and a related decrease in 
population inversion, also leads to decreased peak height. 
This scenario describes the formation of the conductance 
peaks at finite bias. 

Finally we would like to point out the experimental 
significance of the results presented in this section. By 
embedding the quantum dot in e.g. a microwave cav- 
ity the dot can be connected to a photon mode. Due 
to the population inversion this couping might lead to 
stimulated emission of photons, i.e. a maser is formed. 
Previously it has been suggested that such a population 
inversion can be reached in double quantum dot systems, 
where the inversion results from an asymmetric couping 
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of the dots to left and right lead4£ The results of this 
section suggest that an increased bias will modify the 
dot Hamiltonian of two-level quantum dots so that such 
couplings are achieved, as long as the original dot levels 
have different parity. The difference in parity guarantees 
a non-zero dipole element between the dot levels. One 
would expect the dipole element to be larger in two-level 
quantum dots than in double quantum dots due to the 
spatial overlap of the two states. This suggests an in- 
creased observability of population inversion in two-level 
quantum dots. 

VIII. CONCLUDING REMARKS 

A thorough study of the transport through spinless 
two-level quantum dots has been presented. Regimes not 
previously studied in the presence of co-tunneling, such 
as the regimes of finite bias and temperature, as well 
as the transport away from electron-hole symmetry, was 
studied. To present a coherent picture of the transport 
through the system, comparison with previously derived 
results were made. 

Close to degeneracy a suppression of the current was 
found, while a slight detuning of the energy levels re- 
sulted in current peaks. This resulted in a canyon of 
current suppression cutting through both the sequential 
and co-tunneling regime. Different mechanism for the 
suppression were found in the different regimes. In the 
Coulomb blockade region an intuitive explanation to the 
suppression was presented in the low-energy regime, in 
terms of canceling of different tunneling processes due 
to equal occupation of the two levels. In the high-bias 
limit with infinite charging energy a complete suppres- 
sion was found since the electron was trapped in a non- 
conducting state, preventing any transport. These two 
mechanisms for the suppression are independent of the 
lead-dot couplings. When neither of the above mecha- 
nism were present, different parity of the two dot-levels 
had to be assumed to achieve current suppression. For 
such couplings the canyon was found close to degener- 
acy in all four cases (high and low bias, strong and weak 
couplings), independent of the level positions, albeit the 
current did not always entirely drop to zero. 

To quantify the current suppression, simulations were 
performed using the second-order von Neumann tech- 
nique. In the regime where the lead-dot couplings domi- 
nate over temperature and applied bias, i.e. co-tunneling 
processes contribute strongly to the current, a strong sup- 
pression was observed both in the sequential and the co- 
tunneling regime, assuming different parity of the two 



levels. For weaker couplings, where co-tunneling could 
be neglected, only partial suppression was observed in 
agreement with analytical results. 

Away from the electron-hole symmetry point a shift of 
the conductance dip was observed. The shift was shown 
to have three origins: (i) a renormalization of the energy 
levels, (ii) a shift required to give the weakly coupled 
level a higher occupation, and (iii) a shift related to the 
broadening of the levels. 

Finally the conductance peaks were investigated. The 
peaks were found to be intimately related to an increased 
occupation of the upper level. The increased population 
resulted from the bias-induced anticrossing between the 
two orbital levels. Sufficiently high bias resulted in an 
anticrossing large enough to cause population inversion. 
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Appendix A: CURRENT IN THE FIRST-ORDER 
von NEUMANN APPROACH 

In this appendix we will derive expressions for the cur- 
rent assuming that the broadening can be neglected. The 
results are thus of first order in T. In the von Neumann 
formalism the Markov approximation is made within the 
Schrodinger picture, as compared to the Redfield method 
where the Markov limit is taken in the interaction pic- 
ture. This results in slightly different results when the 
levels are non-degenerate4S 

Eq. (|4]) defines t n g(k) as the coupling of a single par- 
ticle state n to state k in lead i. The couplings of the 
corresponding many particle states are given by 

T bae (k)=J2trAk)(b\ll\a), (Al) 

n 

where | a) and | b) are dot states. We use the convention 
that the many particle state b contains one more dot 
electron than a and so on. Assuming energy independent 
couplings, Tb a i(k) = T^, the equation of motion for the 
reduced density matrix is in the first order von Neumann 
method given by4£ 
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ih—w w = (E b - E b ,)w w 
at 



(A2) 



.kt 



■,kl 



Eg' W aa ,T* al J k - J2b" T b"al W b"b'{l - fk 

E k - E v +E a ~ iO+ 

J2b" Tcb"fW b "b'fk ~ J2c> Wcc'T c > b >e(l - f k ) 
' E k -E c + E v + iO+ 



Eg' Tbg'eWg'afk — Eb" w bb"T b "ag{l — fk) rp 



E k -E b + E a + iO+ 



Vat 



\- Eb" Wbb"T* h ,,Jk - Ec' T c>b£ W c>c{l ~ fk) - 

2 , 5 s — — — itti J *i, 



E k -E c + E b - iO+ 



where is the Fermi distribution. 

Going to the continuum limit one replaces the sums 
over k with integrals. The real parts of the integrals con- 
stitute a renormalization of the energy levels. However, 
the results of this Appendix are applied either to electron- 



hole symmetry or to the high/ultra-high-bias case. In 
these situations the renormalization vanishes, and for the 
couplings of Eq. (J6)) the stationary lvN equations read for 
the vacuum state 6 = 0, 6' = 0: 



- w m {f L {E,) + fn(Ei) + a 2 [f L (E 2 ) + f R (E 2 )]} + w„ {2 - f L (E,) - f R (E,)} (A3) 
+ a 2 w 22 {2 - f L (E 2 ) - f R (E 2 )} + aK{w 12 } {f L (E,) + f L (E 2 ) ~ f R {E,) - f R (E 2 )} = 0, 

for the single particle state b — 1, b' = 1: 

™oo {h(Ex) + ME,)} - w n {2 - f L (E,) - /r(£i) + [/l(£ 2 + U) + f R (E 2 + U)]a 2 } (A4) 
+ a 2 w dd {2 - f L (E 2 + U)- f R (E 2 + U)} - aU{w, 2 } {f L (E,) + f L (E 2 + U)- ME,) - f R (E 2 + U)} = 0, 

for the other single particle state b — 2, b' = 2: 

a 2 ^oo Ul(E 2 ) + f R (E 2 )} - w 22 {[2 - f L (E 2 ) - f R (E 2 )]a 2 + f L (E, +U) + f R (E, + U)} (A5) 
+ w dd {2 - f L (E, +U)- f R (E, + U)} - aU{w, 2 } {f L (E 2 ) + f L (E, + U) - f R (E 2 ) - f R {E\ + U)} = 0, 

and for the coherence 6=1, 6' = 2, (6 = 2, V = 1 gives the same information since w 2 i = wl 2 ): 

- {h{Ex) + f L (E 2 ) - ME,) - f R (E 2 )} - ia^ {ME, + U) + f L (E 2 + U) - f R (E, + U) - f R (E 2 + U)} 

- ia^i {ME,) - ME,) + f L (E 2 + U)- f R (E 2 + U)} - ia^f {f L (E 2 ) - f R (E 2 ) + f L (E, +U)- f R (E, + U)} 
- i^[2 - f L (E 2 ) - f R {E 2 ) + f L (E 2 + U) + f R (E 2 + U) (A6) 

+ a 2 (2 - f L (E,) - f R (E,) + f L (E, + U) + f R (E, + U))] + = 0. 

t i 



The equation for 6 = d, b' = d (doubly occupied state) 
does not provide any new information. One must instead 
use the normalization 

w 00 +w,, + w 22 + w dd = 1. (A7) 
At E, = E 2 = — U/2, where wqo = w dd , the equations of 



motion are symmetric in the two levels. This results in 
an equal population of the two states. Furthermore the 
(E, — E 2 )-term of Eq. (|A6j) disappears, implying that the 
coherence is real as remarked in Sec. V. 

Once the system of equations has been solved, the cur- 
rent can be calculated with 40 
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J = y {woolfUEi) + f L (E 2 )a 2 } - w n [l - f L (Ei)] + w n f L (E 2 + U)a 2 - w 22 [l - f L (E 2 )]a 2 + w 22 ! L (E x + U) 
-w dd [l - f L {Ex + U) + [1- f L (E 2 + U)]a 2 ] + aTZ {w 12 } [2 - f L {E x ) - f L (E 2 ) + f L (E x + U) + f L (E 2 + U)}}. (A8) 



We first consider the ultra-high-bias limit, i.e. 
h{E{) = f L (E 2 ) = f L {E x + U) = f L (E 2 + U) = 1 and 
fn(Ei) = f R (E 2 ) = f R (Et +U) = f R (E 2 + U) = 0. In- 
serting these occupations in Eqs. (|A3I) - (|A6I) one obtains 



- u) 00 (l + a ) + wu + a 2 w 22 + 2aK {wi 2 } = 0, (A9) 

w o - w n (l + a 2 ) + a 2 w dd - 2aTZ{w 12 } = 0,(A10) 

a 2 w 00 - w 22 (l + a 2 ) + w dd - 2aK{w 12 } = 0,(A11) 
—\awoo — iawn — iaw 22 — iaw dd 

-i(l + a 2 )u; 12 + -i- — 2 Wl2 = 0.(A12) 
1 1 

The imaginary part of w\ 2 only occurs in the last equa- 
tion and can easily be expressed in terms of the real part 
of W\ 2 . 



X{wi 2 } 



(Et- E 2 )H{w 12 } 



ri(l + a 2 ) 
Eq. (|A12I) can then be rewritten as 



(A13) 



awoo + awn + aw 22 + aw dd 



+ 



1 + a 



2 , (Ei-E 2 ) 



2 1 



Ti(l + a2) 
Solving for the occupations one finds 



Tl{wi 2 }=0. (A14) 



wu = w 22 



4 (i + a 2 ) 2 + (Ex - E 2 )yr( 



, (A15) 



as noted in Sec. VII. Using Eq. (|A8[) the current can be 
calculated as 



J = — [w 00 (l + a 2 ) + w n a 2 
+w 22 + 2aR{wi 2 }] . 

Solving Eqs. (|H|) - (|A12p results in 



(A16) 



It is easy to show that the conductance minimum is found 
at Ei — E 2 where, for a general value of the bias, the 
current is given by Eq. (|2"2")l . 

We next turn to the high-bias limit U = oo, i.e. 
h{Ex) = !l(E 2 ) = 1 and f R (E x ) = f R (E 2 ) = f R {E x + 
U) = f R (E 2 + U) = f L (Ei + U) = f L (E 2 + U) = 0. Here 
it trivially follows from the equation for b = b' = d that 
w dd is zero. Furthermore Eq. (|A3|) is redundant. The 
remaining equations for the occupations are given by 



woo — w n — oJZ {W12} = 0, 
a 2 woo — a 2 w 22 — aTZ {wi 2 } = 0, 
. Wu ■ w 22 
-lawoo - ia— ia ~ 



-l— (1 + a ) + — — w 12 = 0. 

From Eq. (|A20p 1 {w 12 } can be determined 



(A18) 
(A19) 



(A20) 



1{wi 2 ] = - 



2{Ei -E 2 )Tl{w 12 } 



ri(i + <z 2 ) 

Thus in the case U — oo the current is given by 



(A21) 



J = 



hi 



«>oo(l 



£i (l + a^jE^E^/T 2 
h (l + a 2 ) 2 + 3(S 1 -S 2 ) 2 /r 2t 



(A22) 



Once again the conductance minimum is at degeneracy 
but for U — oo there is complete conductance suppres- 
sion independently of the relative coupling strength a, in 
agreement with Refs. [28ll29l 



Appendix B: Renormalization of the dot-levels from 
the Schrieffer- Wolff transformation 



J = 



ri (i 



{E.-Erf/Y 2 } 



2[(l + a 2)2 + ( £ 1 _£ 2 )2 /r 2 ] 



In this appendix we give the derivation of Eq. (|17[) . 
(A17) The renormalization, given by the k = fc'-terms, reads 
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2 (E^^^-^^^+cLc^)) + ^ _^ + ^ l(4c fc , + c\ R C kR ) j dtd! 



/ —at 1 tl\ 



where the factor of 2 accounts for the Hermitean conju- 
gate in Eq. (jlT))) . One sees that the sums are actually 
logarithmically divergent, which of course seems like a 
problem at first sight. However, a common renormal- 
ization of the two levels is just a shift of the zero point 
in the Hamiltonian, that can be neglected. The physi- 
cally interesting quantity is the difference in renormal- 



ization of the two levels. We first assume that level 2 
is fixed, transferring the whole renormalization to level 
1. Replacing the operators with their expectation val- 
ues cfaCk, « {c\ L Ck L ) = fi(E k ) = f(Ek), in the low-bias 
limit, and using Ti(E) = 2nJ2 k t 2 S(Ek — E) and a con- 
stant density of states, the renormalization is given by 



— - t (- p {£>^ (i -'< E *»} -"{£> *- m (B1) 



with V denoting the principal value of the integral. In- 
serting f(E k ) = 1 for E k < and f(E k ) = for E k > 0, 
corresponding to the low temperature limit, and adding 
the constant term a 2 \n\E2/(E2 + U)\ we obtain Eq. (|T7|) 



£i[ln 



E x 



Ei + U 



d\d x 



-a 2 ln 



E 2 




E 2 + U 





(B2) 



This expression can also be derived from first-order 
Redfield which is similar to first order von Neumann, ex- 
cept that the Markov limit is taken in the interaction 
picture, while the Schrodinger picture is used in the von 
Neumann method. Redfield generates an equation of mo- 
tion like Eq. (|A2|) . but with different energy denomina- 
tors in the equations for the off-diagonal terms of the 
reduced density matrix. 



where we have used d\d\ + d\d% — 1, which is applicable 
due to the single occupation of the dot. 
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